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ABSTRACT 


A  modification  to  the  Dantzig  and  Fulkerson  tanker  scheduling  problem 
is  described.  An  insufficient  number  of  vehicles  and  a  utility  associated 
with  each  vehicle  delivery  are  assumed.  The  new  problem  is  shown  to  be 
equivalent  to  a  transshipment  problem,  the  solution  of  which  is  the  same 
as  the  maximal  utility  solution  of  the  modified  tanker  scheduling  problem. 
An  example  is  given. 
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SECTION  I 


INTRODUCTION 

This  paper  treats  a  modification  of  the  tanker  scheduling  problem  solved 
by  Dantzig  and  Fulkerson.  ^ 1  ^  The  problem  was  to  minimize  the  number  of 
vehicles,  M'  ,  needed  to  meet  the  entire  schedule  with  zero  deviation  (i.  e. , 
no  shipment  is  to  be  early  or  late).  Dantzig  and  Fulkerson  were  able  to  show 
that  what  appeared  to  be  a  combinatorial  problem  could  actually  be  formulated 
as  a  Hitchcock  transportation  problem. 

The  problem  considered  in  this  paper  treats  the  case  when  there  are 
insufficient  vehicles,  M  <  M',  *  to  meet  the  entire  schedule,  and,  therefore, 
some  deliveries  must  be  cancelled.  Utilities  are  assumed  for  each  delivery, 
and  the  schedule  is  found  which  maximizes  the  sum  of  the  utilities  of  the 
deliveries  made  with  zero  deviations. 

A  directed  linear  graph,  G'  ,  is  constructed  with  the  property  that  each 
chain  is  a  feasible  vehicle  schedule.  A  set  of  costs  is  determined  by  employ¬ 
ing  a  longest  chain  algorithm.  A  new  problem,  in  the  form  of  the  transship¬ 
ment  problem,  is  solved  to  obtain  the  maximum  flow  solution  with  minimum 
cost.  This  flow  of  value  M  is  shown  to  represent  a  feasible  scheduling  of 
M  vehicles,  which  maximizes  the  total  utility  of  deliveries. 

If  all  utilities  are  positive  and  equal,  then  the  number  of  shipments 
made  with  zero  deviation  is  maximized.  For  simplicity,  it  is  assumed 
that  the  utilities  are  positive  integers. 

(Wherever  positive  integers  are  assumed  in  this  paper,  any  set  of 
positive  rational  numbers  might  also  be  assumed. ) 

♦The  method  can  be  used  also  when  M  =  M1  . 
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SECTION  II 


NOTATION 


Throughout  this  paper,  the  notation  of  Dantzig  and  Fulkerson^  ^  will  be 

adhered  to  wherever  possible.  When  needed,  additional  notation  in  the  spirit 

[2] 

of  Ford  and  Fulkerson  will  be  employed. 

Let  I  =  |  i  |  correspond  to  the  set  of  ports  at  which  shipments  can 
originate;  let  J  =  J  j  j  correspond  to  the  set  of  ports  at  which  a  delivery 
is  to  be  made.  *  Associated  with  the  sets  I  and  J  are  sets  of  positive 
integers  A  =  j  a^  j  ,  the  transit  time  from  i  to  j  (including  loading  and 
unloading),  and  B  =  j  b^  j  ,  the  return  travel  time  from  ]  to  i  . 

The  kth  shipment  from  i  to  j  requires  m^  >  0  vehicles,  leaving 

k 

(more  accurately,  starting  to  load  at )  port  i  at  time  t^  and  arriving 
(more  accurately,  completely  unloaded)  at  port  j  at  time 


T*4 

1J 


‘ij  +  atJ 


(1) 


k  k  k  k  k 

Each  partial  shipment  has  a  utility,  d^  v^  s  m^  v^,  where  v^  is 

k 

the  utility  of  a  single  vehicle  delivery  and  d  are  the  number  of  vehicle 

k  k  +  1 

loads  delivered.  It  is  assumed  that  t^  <  t^ 

A  directed  graph  G'  ,  is  defined  consisting  of  sets  of  nodes,  s'  ,  s  , 

X,  Y,  t.  Each  node  xcXis  denoted  by  a  unique  pair  of  positive  integers 

k 

(of ,  i).  A  node  x  ~  (a  ,  i)  exists  if,  and  only  if,  at  least  one  t  =  a.  In 
like  manner,  each  node  y  e  Y  is  denoted  by  a  unique  pair  of  positive  integers 
(fi  ,  J).  A  node  y  ~  (fi  ,  j)  exists  if,  and  only  if,  at  least  one 


*Note  that  the  intersection  of  sets  I  and  J  is  not  necessarily  null. 
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R  R 

Tij  =  *ij  +  aij  =  fi  •  Further,  a  source  node,  s'  ,  corresponding  to  an 
artificial  source  of  vehicles,  a  sink  node,  t ,  and  a  node  s  are  defined. 


For  each  pair  of  nodes,  x  ~  (a  ,  i)  and  y  ~  (J3  ,  j),  a  set  of  P 

k  k  +  1 

parallel  arcs,  (x  ,  y)  ,  is  defined  if,  and  only  if,  t^  =  =  .  .  . 

=  t^j+  P  ”  1  =a  •  (N°te  also,  that  =  T^+  1  =  .  .  .  =  T^+  P  1  =/?, 

R  +  p  R  +  p 

since  =  t^  +  a^  . )  The  capacity  and  utility  of  each  arc  are 

defined  as  c  (x  ,  y)  =  +  P  ”  1  and  v  (x  ,  y)  =  +  p  "  1  respectively, 

P  ij  P  ij 

P  =  1.  2 . P  . 


Further,  for  each  pair  of  nodes,  y  ~  (fi  ,  j)  and  x  ~  (o; ,  i),  an  arc, 

(y  ,  x)^  ,  is  defined  if,  and  only  if,  a  2  M  ;  with  capacity,  c^  (y  ,  x) 
=  «° ,  and  utility,  v^  (y  ,  x)  =  0.  Finally,  arcs,  (s  ,  x)^  ,  and  arcs, 

(y  ,  t)1>  are  defined  for  x  e  X  and  y  f  Y,  with  c^  (s.  x)  =  (y  ,  t) 

=  «  and  Vj  (s  ,  x)  =  Vj^  (y  ,  t)  =  0. 
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section  m 


EXAMPLE 


Consider  a  tanker  scheduling  problem  involving  two  origins  and  two 
destinations.  The  shipment  departure  times  are  given  in  a  rectangular 
array,  one  row  for  each  origin,  and  one  column  for  each  destination. 

In  each  space,  (i  ,  j)  ,  is  the  sequence  of  shipment  departure  times: 


j  =  1 


i  =  1 
i  =  2 


1  =  2 


1 

9  ,  9 

9 

1.  5 

kij 


The  number  of  vehicle  deliveries  comprising  each  shipment  is  given 
below: 

j  -  1  j  *  2 


i  =  1 
i  =  2 


1  ,  1 


1  ,  1 


=  <  m. 


The  utility  of  each  vehicle  delivery  is: 

J  =  1  j  =  2 

I 

i  =  1 


i  =  2 


2  ,  1 


2  ,  1 


=  <v 


ij 


For  this  example,  assume  that  the  transit  times  are: 
1=1  1=2 


i  =  1 
i  =  2 


*U 


=  b, 


"ij 
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From  Equation  (1)  we  construct  the  table  of  shipment  delivery  times 
and,  thus,  the  time  when  vehicles  are  available  for  reassignment: 

j  =  1  3  =  2 

i  *  1 

i  =  2 

Consider  the  shipment  from  port  1  =  1,  departing  at  time  t^  =  1, 
consisting  of  m^  =  2  vehicle  loads,  and  arriving  at  port  j  =  1  at 

tJj  =  8.  Each  load  has  utility  v^  =  2,  hence  the  total  utility  of  this 
shipment  is  4.  For  this  delivery,  construct  the  nodes,  x  ~  (a  *  1  ,  i  =  1) 

and  y  ~  (/3  =  8,  j  =  1),  and  the  arc,  (x  ,  y)  ,  with  capacity,  c  (x  ,  y) 

1  1  1 
=  m^  =  2,  and  utility,  v  (x  ,  y)  =  v  ■  2.  In  a  similar  manner,  the 

arcs  and  nodes  corresponding  to  the  other  shipments  are  constructed  as 

shown  in  Figure  1.  (Note  the  two  parallel  arcs  between  x  ~  (a  =  9  ,  i  =  1) 

and  y  -  (/3  ■  12,  j  =  2)  ,  which  correspond  to  two  vehicle  loads  having 

different  utilities. ) 

Arcs  are  now  added  to  correspond  to  feasible  reassignments  of  vehicles 
from  a  destination  port  to  an  originating  port.  For  example,  consider  the 
node,  y  ~  (/?  =  3,  j  =  2).  If  this  shipment  is  made,  then  the  vehicle 
could  be  reassigned  to  any  origin,  x  ~  (a  ,  i)  ,  such  that  a  z  / 3  +  , 

specifically,  to  either  y  ~  (a  =  5,  i  =  2),  (a  =  9,  i  =  2)  ,  or 
(a  =  9  ,  i  =  1).  No  other  origins  are  feasible  with  respect  to  transit 
times.  Each  of  these  arcs  has  infinite  capacity  (recall  arcs  terminating 
at  y  ~  (/J  ,  j)  have  finite  capacity)  and  zero  utility.  Figure  2  illustrates 
the  construction  so  far. 

Since  any  shipment  may  be  made  by  aBBigning  a  new  vehicle  instead 
of  reassigning  an  existing  vehicle,  arcs  are  added  from  s  to  each  node 
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FIGURE 


FIGURE  2 

ADDITION  OF  REASSIGNMENT  ARCS 


x  ~  (a  ,  i).  Each  such  arc  has  infinite  capacity  and  zero  utility.  In  a  similar 
manner,  we  add  arcs  of  infinite  capacity  and  zero  utility  from  each  node, 
y  ~  (0  ,  j),  to  t  which  correspond  to  not  reassigning  the  vehicle  after  a 
delivery  is  made.  Finally,  we  add  node,  s'  ,  and  arc,  (s'  ,  b)^  ,  with 
capacity  equal  to  the  number  of  vehicles  available,  M  ,  and  zero  utility. 

The  completed  graph  is  illustrated  in  Figure  3. 

From  the  method  of  construction,  it  is  seen  that  there  is  a  one-to-one 
correspondence  between  shipping  schedules  for  M  vehicles  and  flows  of 
value  M  from  s'  to  t  in  the  constructed  graph,  G'  . 

G'  is  an  acyclic  graph,  since  arcs  from  x  ~  (a,i)  to  y  ~  (/?  ,  j)  have 
the  property  that  /3  =*  a  and  arcs  from  y  ~  (P  ,  j)  to  x  ~  (cv,  i)  have  the 
property  that  a  >  1 3  ,  and ,  hence,  for  any  chain,  the  times,  a  or  p  ,  are 
monotonically  increasing  and  no  cycle  can  occur. 
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COMPLETED  GRAPH 


SECTION  IV 


DETERMINATION  OF  ARC  COSTS  AND  PROBLEM  SOLUTION 

Since  it  has  been  assumed  that  M  ,  the  number  of  vehicles  available, 
is  less  than  or  equal  to  M'  ,  the  minimum  number  of  vehicles  required  to 
meet  all  shipping  dates,  a  method  to  find  a  flow  of  value  M  which  maximizes 
the  total  utility  of  all  shipments  delivered,  i.  e. ,  use  all  available  vehicles, 
is  desired. 

Next,  a  method  is  presented  of  assigning  positive  arc  costs  such  that 
the  minimum  cost  flow  of  value  M  corresponds  to  the  maximum  utility 
solution  to  the  original  problem.  * 

A  simple  modification  of  the  shortest  chain  algorithm  can  be  made 

which  will  find  the  longest  chain  from  any  node  to  a  specific  node  in  an 
r  3i  r4i 

acyclic  graph.  J  1  Denote  ?r(e)  as  the  length  of  the  longest  chain 
from  any  node,  e  ,  in  G'  to  node  t ,  where  the  length  of  an  arc  is  v  (e  ,  f). 
At  the  termination  of  the  longest  chain  algorithm,  the  ir  (e)  correspond  to  the 
node  numbers. 

Define  a  (e  ,  f)  ■  ir  (e)  -  n  (f)  -  v  (e  ,  f)  as  the  cost  associated 
P  P 

with  all  arcs  (e  ,  f)  ,  of  G'  .  Figure  4  shows  the  costs  and  node  numbers 
associated  with  the  graph  G'  .  Each  node  has  been  assigned  an  identifying 
number  for  future  reference. 


♦Since  G'  is  acyclic,  an  alternative  solution  method  ib  to  assign  arc  costs 
equal  to  the  negative  of  the  utility  of  the  shipments  over  each  arc  and  to  find 
a  minimum  cost  flow  in  G'  .  One  could  use  the  "Out-of  Kilter”  algorithm 
described  in  (2)  or  the  primal  method  described  in  (5). 
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(C=  2t  V=2 ) 
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FIGURE  4 

DETERMINATION  OF  ARC  COSTS 


THEOREM 


The  total  cost  associated  with  a  unit  chain  flow,  in  the  q  chain, 

[  s'  ,  (s'  ,  s)  ,  s  ,  (s  ,  x1)i  ,  x1  ,  (xL  ,  y1)p  ,  yx,  .  .  .  ,  \  .  yn)p  . 

y  .  (y  ,  t),  ,  t]  ,  in  the  directed  graph  G'  ,  is  given  by: 
n  n  1 


C  =  T(S')  - 

q 


(Xq  '  V 


v  (x  ,  y  ) 
Pq  q  q 


for  all  nodes  x  c  X  ,  y  «  Y  in  the  chain. 

q  q 


PROOF 


The  total  chain  cost  is: 


Cq  =  “l  ,S'  •  S)  "  :'l  (S  '  V  ^  aPl  (X1  •  yi>  +  ai  (yi  *  V 


-  .  .  .  -  a  (y  .  t) 
1  "  « 


C  =  T(S')  -  7T(S)  -  V  (S'  ,  S)  +  T(S)  -  7T(X  )  -  V  (S  ,  X  ) 

q  1  111 

+  r  (Xj )  -  ”  ty^)  -  vp  (Xx  ,  y2)  +  xiyj  -  *(x2) 


-  v!  <yx  •  V  +  •  •  •  +  ■  ’r(t>  "  V1  (yn  ’ 


but,  the  only  nonvanishing  utilities  are  v  (x  ,  y  )  and  n(t)  =  0. 


Hence,  letting 


pm  m  m 


be  the  summation  over  all  arcs  in  the  q  chain 


(xq  '  yq> 
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C  =  IT  (S1)  - 

q 


Q.e.d, 


7  Vn  <x„  •  y„> 
u  p_  q  q 

<x„  •  yj 

q  q 


Note  in  the  theorem,  the  total  cost  of  a  unit  flow  in  chain  q  is  r  =  ff(s') 

q 

-  (the  utility  of  loads  delivered  over  chain  q).  The  total  cost  of  a  chain  flow 
of  value  M  is  the  sum  of  the  costs  of  the  M  unit  chain  flows,  and  is  given  by: 
M  M 

C  = 


l 


C  =*  Mjt 

q 


<8')-  I  I  Vv 


V 


q  =  1 


q  =  1  <vV 


C  =  M  it  (s')  -  Y,  Z  dfj  vfj  =  Mff(s,)  "  V 


i  j  k 


where, 


i  j  k 


For  M  fixed,  M7r(s')  is  a  constant  and  the  flow  of  value  M  which  mini¬ 
mizes  the  total  cost  also  maximizes  the  total  utility  of  the  loads  delivered  by 
M  vehicles. 


The  transshipment  problem  (general  minimum  cost  flow  problem)  is 
to  find  a  maximum  flow  solution  having  the  minimum  cost.  Since  the  maxi¬ 
mum  flow  for  our  problem  is  M  =  c}  (s'  ,  s)  ^  M'  ,  the  transshipment 
problem  is  to  find  a  maximum  flow  of  value  M  having  minimum  cost,  and 
therefore  maximum  utility.  * 


* 

If  it  were  not  for  c^  (s'  ,  s)  ,  a  flow  of  value  M'  or  greater  would  be 
feasible. 
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The  above  construction  of  the  graph,  G'  ,  has  transformed  the  tanker 

scheduling  problem  with  a  limited  number  of  vehicles  into  a  transshipment 

problem.  Most  transshipment  algorithms  produce  the  optimum  flow  in  node- 

arc  form;  however,  since  G'  is  an  acyclic  graph,  the  optimum  flow  can 

[2] 

easily  be  transformed  into  arc-chain  form. 

METHOD 

(1) 

(2) 

(3) 

(4) 

(5) 


NUMERICAL  RESULTS  FOR  THE  PREVIOUS  EXAMPLE 

For  M  =  1  ,  the  transshipment  algorithm  terminates  with  unit  flow 
in  the  following  arcs:  (s',  s)1  ,  (s  ,  x^  ,  (xg  ,  y^  ,  (yg  ,  x^  , 

(x3  »  '  ^3  •  X4>1  •  <x4  •  ^4*!  •  ^4  •  •  1116  minimum  cost  is  zero 

and  the  utility  of  deliveries  made  is  v(s')  -0  =  5.  Transformation  of  this 
flow  yields  the  following  schedule  for  the  single  vehicle  (which  is  unique). 


OF  SOLUTION 

Construct  the  directed  graph,  G'. 

Use  a  longest  chain  algorithm  to  find  w  (e)  ,  for  all 
nodes,  e  , 

Calculate  the  arc  costs,  a  (e  ,  f)  =  ir(e)  - 
*  (f)  -  vp  (e  ,  f)  . 

Solve  this  transshipment  problem  using  an  appropriate 
algorithm  to  obtain  the  maximum  utility  solution. 

Construct  individual  vehicle  schedules  from  the  optimum 
arc  flows  in  G'  by  transforming  them  into  M  unit  chain 
flows,  each  corresponding  to  an  optimum  vehicle  schedule. 
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ft 

1 


X 


Leave 

From 

Arrive 

At 

Utility 

a  =  1 

i  =  2 

TS> 

II 

CO 

j  =  2 

2 

a  =  5 

i  =  2 

3=7 

J  =  2 

1 

a  =  9 

i  ^  2 

o 

f-H 

II 

j  =  1 

2 

V  =  5 

The  solution  for 

M  =  2 

vehicles  is: 

Vehicle  1 

Leave 

From 

Arrive 

At 

Utility 

o  =  1 

i  =  1 

CO 

II 

*33 

J  =  1 

2 

o’  =  9 

i  =  2 

O 

II 

j  =  1 

2 

II 

> 

Vehicle  2 

Leave 

From 

Arrive 

At 

Utility 

o  1 

i  =  2 

J  =  :t 

j  =  2 

2 

a  =  9 

i  =  1 

if  =  12 

j  =  2 

2 

For  two  vehicles,  the  optimum  schedule  is  unique  and  does 
not  include  the  optimum  schedule  for  the  case  when  M  =  1,  which  serves 
to  show  that  a  universally  optimum  schedule  does  not  exist. 

All  deliveries  can  be  made  when  M  =  M'  =  5.  If  M  =  6,  the 
algorithm  would  use  all  six  vehicles  and  should  not  be  used  for  values  of 
M  >  M'  .  M'  can  be  determined  either  from  the  Dantzig  and  Fulkerson  ^  1  ^ 

algorithm  or  by  using  this  algorithm  iteratively,  M  =  1  ,  2 . M' , 

and  determining  M'  as  the  minimum  value  M  for  which  all  shipments 
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i  J  k 

would  not  be  done  unless  a  complete  parametric  analysis  is  desired. 

An  alternative  procedure,  employing  the  "Out-of- Kilter"  algorithm, 
would  assign  a  small  positive  cost  v^s',  s)  =  e  >  0  to  the  arc  (s',  s^  ; 
assign  costs  v'  (x  ,  y  )  =  -  v  <x  ,  y  )  to  the  delivery  arcs;  and 

Pq  q  q  Pq  q  q 

determine  the  minimum  cost  flow  solution  in  the  graph  G' .  This  proce¬ 
dure  may  be  employed  even  when  M  =»  M1 . 
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